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ABSTRACT 

Pain is more than an unpleasant sensory experience 
associated with actual or potential tissue damage: it 
is the most common reason for physician consult- 
ation and often dramatically affects quality of life. 
The management of pain is often difficult and new 
targets are required for more effective and specific 
treatment. SuperPain (http://bioinformatics.charite. 
de/superpain/) is freely available database for 
pain-stimulating and pain-relieving compounds, 
which bind or potentially bind to ion channels that 
are involved in the transmission of pain signals to 
the central nervous system, such as TRPV1, TRPM8, 
TRPA1, TREK1, TRESK, hERG, ASIC, P2X and 
voltage-gated sodium channels. The database 
consists of ~8700 ligands, which are characterized 
by experimentally measured binding affinities. 
Additionally, 100000 putative ligands are included. 
Moreover, the database provides 3D structures of 
receptors and predicted ligand-binding poses. 
These binding poses and a structural classification 
scheme provide hints for the design of new anal- 
gesic compounds. A user-friendly graphical inter- 
face allows similarity searching, visualization of 
ligands docked into the receptor, etc. 

INTRODUCTION 

Ion channels are proteins forming a pore that allows the 
flow of ions across membranes. Ion channels are voltage 
or ligand gated. Some of these proteins help nerve cells to 
transmit pain signals to the central nervous system and are 
therefore promising targets for the development of pain 
therapeutics. 

Transient receptor potential channels (TRPs) are a 
family of 28 human cellular ion channels, varying in 



homology to each other but all with six transmembrane 
regions in common, and are nonselectively permeable to 
cations. There are seven subfamilies that can be divided 
into two groups. Group 1 includes TRPC, TRPV, TRPA, 
TRPM, TRPN, and group 2 comprises TRPP and 
TRPML (1). Compounds perceived as hot stimulate the 
vanilloid receptor (TRPV). At the same time, this receptor 
plays a crucial role in pain mediation and is therefore an 
interesting drug target. Known pungent chemicals with 
high receptor affinity such as capsaicin or resiniferatoxin 
were lead structures in drug development toward desensi- 
tization (2). Recent research focuses on novel analgesic 
mechanisms like positive allosteric modulation (3) and 
broad-spectrum TRP antagonists (4). 

The mechanism of cold-induced analgesia was unclear 
until the discovery of TRPM8. This receptor is also known 
as the cold or menthol receptor (5). It is activated by cold 
temperatures and cooling agents, such as menthol or icilin, 
allowing the entry of Na + and Ca 2+ to the cell (6). Two 
modulating mechanisms are generally discussed. Whereas 
antagonists physically block the receptor for cold and 
menthol, agonists activate TRPM8 and generate a 
cooling sensation. Selective ligands could be used as a 
new generation of analgesic drugs in neuropathic pain 
(7,8). 

TRPA1 plays a key role in chemical sensing in the in- 
flammatory pain pathway. Many small molecules, 
including ingredients like wasabi, horseradish, garlic and 
mustard oil, can activate the channel. Recently, it has been 
shown that desensitizing TRPA1 could help in the treat- 
ment of neuropathic pain (9,10). 

The human Ether-a-go-go Related Gene (hERG) 
channel or KCNH2 is a voltage-gated potassium 
channel. It has been the focus of pharmaceutical 
research for years because the inhibition of hERG potas- 
sium channels by drugs can lead to cardiac arrhythmia 
(11). In 2010, Stary and colleagues published a 
homology model (12) and a group from Canada 
investigated structural mechanisms of state-dependent 
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drug binding (13). The different binding sites for ions and 
compounds are defined and the International Conference 
of Harmonization set up guidelines for drug development. 
Therefore, new screening methods for the prediction 
of drug liability to the hERG channel have been 
developed (14). 

Voltage-gated sodium channels are activated through 
action potential firing. They represent the target for 
local anesthetic agents. It is challenging to find selective 
inhibitors of sodium channels in the pain pathway. 
Genotyping of families suffering from congenital indiffer- 
ence to pain identified mutations in gene coding for 
Navl.7 channels (15). A lot of research is performed 
to find specific blockers to the treatment of pain and 
epilepsy (16). 

There are some complementary resources on ion 
channels. Multiple analysis for voltage-gated potassium 
channels from different species is compiled in the 
voltage-gated potassium channel database VKCDB (17). 
Protein-protein interactions (PPIs) are the focus of the 
TRIP database (18): Shin and colleagues manually 
curated 653 PPIs for mammalian TRP channels. 
MoleOnline 2.0 (19) is a web server that provides inter- 
active channel analysis to identify active sites. 

The International Union of Basic and Clinical 
Pharmacology (IUPHAR) database provides information 
on human and rodent receptors (20). 

There is a need for a specific resource for pain-relieving 
compounds targeting ion channels and their 3D homology 
models so researchers can identify targets and putative 
ligands. 

MATERIALS AND METHODS 

Here, we shortly describe our methods used. If you are 
interested in further details, please read the 'About' page 
on our Web site. 

Compounds 

To identify pain-related receptors and ligands, scientific 
literature was screened via text mining and subsequent 
manual evaluation. Therefore, we first downloaded 
Medline/PubMed data from the NCBI FTP site in xml- 
format. Using the search engine library Apache Lucene 
(http://lucene.apache.org) and a tool kit for processing 
text with computational linguistics (http://alias-i.com/ 
lingpipe), the data was indexed. The search engine, 
written in Java, dynamically queries the indexed data 
and results in a structured query language (SQL) file con- 
taining the textmining hits. An example for a query for the 
literature search is as follows: [lidocaine (TI) AND ic50 
(TI)] OR [lidocaine (abstract) AND ic50 (abstract)]. 
Found hits were put in a preliminary database and are 
displayed with keywords highlighted in different colors. 
In the manual evaluation process, confirmed experimen- 
tally determined affinities were moved to the final 
database. In some cases, the full text of the article was 
checked. Further information regarding the molecules 
was retrieved from the PubChem database (21). 
PubChem is a freely available database, which is 



provided by the National Center for Biotechnology 
Information (NCBI). It contains detailed information on 
~30 million compounds. All those compounds and add- 
itional information were put in a MySQL database. 
Additional affinities were found via a search in 
BindingDB (22). BindingDB is a database on experimen- 
tally determined protein-ligand interactions. It provides 
~782 000 affinities for 6500 protein targets. 

Compound clustering 

The experimentally determined ligands had to be clustered 
regarding their similarity to each other by means of a K- 
means algorithm. As in the K-means algorithm, the 
number of clusters has to be defined the algorithm was 
slightly modified. To ensure that the most similar com- 
pounds are members of one cluster, a neighbor-joining 
algorithm was used. The R package with heatmap.2 was 
used to display the compound similarities in a heatmap. 

In silico screening 

Putative ligands were indentified through a similarity 
search. Therefore, the structural fingerprints of each 
ligand were calculated and used for a similarity search 
within the PubChem Compound database. The similarity 
between experimentally determined ligands and potential 
ligands was calculated with the Tanimoto coefficient. 
Depending on the chemical topological properties it 
gives values between zero (no similarity) and one (identi- 
cal). Compounds with a coefficient >0.90 were classified 
as putative ligands and included in the database. 

Receptor structure 

Currently, no radiographic crystallographic structures are 
available for most of the receptors, but there is some in- 
formation on the active sites of some of the molecules. To 
perform dockings a homology model for TRPV1 had to 
be created. The template structure for the homology 
modeling of TRPV1 and the models of P2X and ASIC 
were retrieved from the Protein Data Bank (PDB). The 
PDB is a freely available database on 3D structures of 
proteins und nucleic acids (23,24). Models are download- 
able and provide coordinates of each atom within the 
molecule. A homology model of the hERG channel was 
retrieved from a research group in Germany (12). 

Docking 

To obtain an understanding of binding mechanisms, the 
homology model of TRPV1 and the models of P2X, ASIC 
and hERG were imported into Accelrys Discovery Studio 
(25). About 6000 ligands were docked into the binding 
sites with the integrated Docking Tool LibDock. 
LibDock is a high-throughput docking algorithm. Based 
on polar and apolar interaction sites up to 1000 ligand 
conformations were positioned in the binding site. The 
five best-ranked poses regarding energetic conditions are 
displayed on the Web site. 'JSmol' was used to implement 
a molecular viewer, which is JavaScript based. 
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Web site 

The Web site is based on PHP (http://www.php.net/); web 
access is enabled through Apache HTTP Server (http:// 
httpd.apache.org/). We recommend a recent version of 
Mozilla Firefox or Google Chrome; alternative browsers 
like Microsoft Internet Explorer and Apple Safari were 
tested with some configurations. JavaScript must be 
enabled, as it is required for all features of the site. 

RESULTS 

Compound database 

A total of 100000 putative channel modulators were 
included of which 12 000 are purchasable and 8700 
ligands have measured binding affinities. 

Compound clustering 

A structural clustering of the 8700 compounds with an 
internal similarity (Tanimoto) above 0.7 was performed 
and resulted in clusters with member sizes up to 52 com- 
pounds. For a better visualization, the member size was 
limited to 30. The neighbor-joining algorithm resulted in 
684 clusters. The similarities are displayed in interactive 
heatmaps. 

Receptor structure 

The homology model of TRPV1 was based on the known 
structure of the Kvl.2 potassium channel, which exhibits 
significant similarity and a high resolution (2.4 •). The 
structure was downloaded in PDB format (PDB-ID: 
2R9R). Homology modeling was carried out in accord- 
ance with Fernandez-Ballester et al. (26) using the 
SWISS-MODEL server alignment mode. PyMOL was 



used to create a tetramer and manual refinement was 
carried out with Accelrys Discovery Studio. 

Docking 

The docking of 6000 TRPV1, P2X, ASIC and hERG 
ligands in > 1000 different binding poses revealed five 
binding poses for each ligand. Figure 1 shows capsaicin 
in the binding site of TRPV1 in an active interactive view. 
These binding poses were chosen regarding energetic con- 
ditions. The results are embedded into the database. Each 
ligand can be displayed in different manners, as well as the 
receptor itself (ball and stick, spacefill, etc.). Holding the 
left mouse button while moving the mouse leads to a 
rotation. A right click brings up a drop-down menu for 
zooming, measurements, color adjustments, etc. Safari or 
Tablet users can use their common gestures to navigate 
through the docking results. Each binding pose can be 
exported as a graphic (PNG) or file (MOL). 

Usage of the database 

There are different ways to use the database and browse 
through the data. Figure 2 summarizes the main 
functionalities of the Web site. 

'Receptor' holds information on pain-related ion 
channels. The homology model of TRPV1 and docking 
results can be found there. 

There are two search boxes in the 'Compounds' section. 
Specific compounds can be found through a property 
search by typing in the name, PubChem ID, IUPAC or 
SMILES. It is also possible to search for a specific target 
such as TRPV1 or hERG or to select certain features like 
purchasability, IC50, EC50, molweight, rotatable bonds, 
etc. On the results page, information on compounds 
including vendors is listed in a table. This table can be 
used as a starting point for a similarity search or to 
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Figure 1. Capsaicin in TRPV1. Capsaicin was put into the homology model of TRPV1 using a high-throughput docking algorithm. Each docked 
ligand is shown in the interactive view. 
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Figure 2. Main functionalities of 'SuperPain'. The database provides 
information on pain-related ion channels, as well as 3D structures 
and interactive docking results. The 'Compound' search enables users 
to find ligands regarding properties or features (name, SMILES, target, 
IC50, molweight, purchasability, etc.). The 'Results' page allows to 
perform similarity searches or to view docking results or to see 
heatmaps with similar structures to compare Tanimoto coefficient or 
affinities. 



compare structures directly by viewing the compound 
clusters. In the 'Cluster' section all clusters are listed 
with one compound of the cluster. A search box enables 
quick access to a distinct cluster. Clicking on a compound 



also enables browsing through different clusters. The 
heatmaps are diagonally divided to compare similarity 
(Tanimoto) and affinity (IC50). Mouseover shows the 
structure of the compared compounds. Clicking leads to 
a comparison page with information about (co-) affinities, 
as well as chemical information and vendors. 



DISCUSSION 

SuperPain aims at providing a comprehensive resource on 
ligands for pain-related ion channels. 

Acid sensing ion channels (ASICs) belong to the 
degenerin-epithelial sodium channel superfamily and are 
voltage independent. They are activated by extracellular 
acidosis and are involved in different processes, such as 
taste, mechanosensation and nociception (27). Recently, it 
has been shown that mambalgin, a three-finger peptide 
from the venom of the black mamba, suppresses pain in 
mice without toxicity and fewer side effects than morphine 
(28). These findings show that ASICs are promising 
targets in the treatment of pain, especially to avoid addic- 
tion problems (29). 

In traditional Chinese medicine, tetramethylprazine, 
sodium ferrulate and puerarin are used in the therapy of 
pain. It has been found that these compounds target some 
P2X receptors. These ligand-gated ion channels open in 
response to binding of adenosine 5'-triphosphaste (ATP). 
Advances in radiographic crystallography allow in silico 
ligand docking and P2X receptors are becoming thera- 
peutically important drug targets (30). 

Further potassium channels of interest are the TWIK- 
related spinal cord potassium channel (TRESK) and 
TREK-1. TRESK is a two-pore domain potassium 
channel that is mainly expressed in dorsal root and 
trigeminal ganglia (31). This channel is involved in acute 
and chronic pain and plays an important role in migraine 
pathogenesis, which makes it a promising target (32). 
TREK-1 is an interesting target because of its co-expres- 
sion with TRPV1 and its involvement in polymodal pain 
perception (33). The present intensive research on 3D 
structures will improve the development of channel modu- 
lators. For example, there are different binding sites for 
agonists and antagonists (34). The binding sites of the ion 
channels are rather large, which allows binding of a 
variety of compounds if certain physiochemical property 
conditions are fulfilled. This is reflected by the great diver- 
sity of the pain-related compounds. They show a low 
mean similarity of 0.4 compared with other targets like 
PARP (0.6). It is also important to find out more about 
the interactions or co-inhibition of different receptors for 
pain. In the database, they can be found by choosing 
'MultiTarget' from the 'Compound Properties' search 
box. For example, menthol is known to be an inhibitor 
of TRPM8 with a half maximal effective concentration 
(EC 50 ) of 29 000 nM. At the same time, it has also inhibi- 
tory potency on TRPA1 with an EC 50 of 28400nM. 
Cannabigerol is a compound that occurs naturally in 
hemp strains. In contrast to other cannabinoids, it is 
nonpsychoactive because it is not only a a 2 -adrenergic 
receptor agonist, but also a mild CB, receptor antagonist 
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(35). Although the compound is not well-studied, it has 
been found to lower the intraocular pressure (36) and 
there is some research on treating inflammatory bowel 
diseases with cannabigerol (37). 'SuperPain' stores 
affinities for TRPA1, TRPM8 and TRPV1 for the 
compound. These experimentally determined co-affinities 
suggest that these multi-target compounds might be 
promising for the development of pain therapeutics (38). 

Availability 

'SuperPain' is publicly available via http://bioinformatics. 
charite.de/superpain and should be used under the terms 
of the Creative Commons Attribution-Noncommercial- 
Share Alike 3.0 License. 
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